Neutron and ARPES Constraints on the Couplings of the 
Multiorbital Hubbard Model for the Pnictides 
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The results of neutron scattering and angle-resolved photoemission experiments for the Fe-pnictide parent 
compounds, and their metallic nature, are shown to impose severe constraints on the range of values that can 
be considered "realistic" for the intraorbital Hubbard repulsion U and Hund coupling / in multiorbital Hubbard 
models treated in the mean-field approximation. Phase diagrams for three- and five-orbital models are here 
provided, and the physically realistic regime of couplings is highlighted, to guide future theoretical work into 
the proper region of parameters of Hubbard models. In addition, using the random phase approximation, the 
pairing tendencies in these realistic coupling regions are investigated. It is shown that the dominant spin-singlet 
pairing channels in these coupling regimes correspond to nodal superconductivity, with strong competition 
between several states that belong to different irreducible representations. This is compatible with experimental 
bulk measurements that have reported the existence of nodes in several Fe-pnictide compounds. 

PACS numbers: 74.70.Xa, 74.20.-z, 74.20.Rp 



I. INTRODUCTION 

An exciting new area of research has recently opened with 
the discovery of superconductivity in the Fe pnictidesji^ 
These materials present many similarities with the high tem- 
perature superconductors based on copper, but also several 
differences. Both families have layered structures suggesting 
that theoretical studies in two dimensional lattices should be 
able to capture the essence of the pairing mechanism. In the 
undoped limit, both types of compounds are magnetic, with 
wavevector {n, 0) in the case of the pnictides**'^ and (tt, n) for 
the cuprates, in the notation of the square lattice defined by 
Fe or Cu. However, a crucial difference is that the undoped 
parent compound is an insulator for the cuprates, while it is a 
(bad) metal for the pnictides. This fact already suggests that 
the regime of a large Hubbard coupling U, widely used in the 
context of the cuprates, may not be appropriate for a theoreti- 
cal description of the pnictides. 

Adding to the complexity of this problem, it is also clear 
that a theoretical study of pnictides cannot rely on just one 
orbital, as in the case of the cuprates, but it needs a multior- 
bital approach.'" In fact, to study the magnetic and supercon- 
ducting properties of the pnictides, a considerable effort has 
already started using multiorbital model Hamiltoniansjii"— 
Having to consider multiple orbitals severely restricts the 
available tools to carry out unbiased computer-based inves- 
tigations of Hubbard multiorbital models. As a consequence, 
several studies have been restricted to mean-field approxima- 
tions. Fortunately, this is not a drastic limitation for the case of 
the undoped systems, since similar mean-field approximations 
for the cuprates are known to capture qualitatively the essence 
of the magnetic states,^^ In fact, recent mean-field-based ef- 
forts for the pnictides have already reported the presence of 



a state that is simultaneously metallic and magnetic.'^ How- 
ever, the study of the intermediate coupling regime, where the 
magnetic-metallic state was found, establishes a considerable 
challenge to theory since an intermediate region of couplings 
is often more difficult to analyze than either extreme of large 
or small U. In addition, a multiorbital Hubbard approach 
needs at least two couplings: the on-site intra-orbital Hubbard 
repulsion U and the on-site Hund coupling J. (A third pa- 
rameter, the on-site inter-orbital repulsion U' , is then defined 
by the well-known relation U - U' + 2J arising from sym- 
metries in orbital space. ■'^) Having two couplings increases 
further the complexity of the analysis and the comparison be- 
tween different approaches since there is at present no univer- 
sally accepted range of U and J that is considered realistic by 
the community of experts. In fact, it would be quite desirable 
to restrict the values of the U and J couplings used in the liter- 
ature to a much narrower range, where qualitative agreement 
with experiments is observed. 

In this publication, our goal is to use experimental neutron 
scattering and photoemission data for the undoped pnictide 
parent compounds, supplemented by their well-known metal- 
lic properties, to establish lower and upper bounds on the cou- 
plings U and J of multiorbital Hubbard models. By focusing 
on a more restricted set of couplings, our results will guide fu- 
ture theoretical efforts into a realistic regime for the pnictides. 
Our calculations are based on the previously used and tested 
mean-field approximation,'** allowing us to calculate a variety 
of observables that are then compared against experimental 
results to establish the proper ranges for U and J. 

To caiTy out this theory-experiment comparison the focus 
here is on the results of two powerful experimental techniques. 
One of them is neutron scattering. In these experiments, a 
magnetic order with wavevector (n, 0) (in the Fe square-lattice 



notation) was observed at low temperatures.^'^ This particular 
wavevector can be theoretically accommodated by merely us- 
ing proper tight-binding model Hamiltonians that reproduce 
the well-known paramagnetic Fermi surface (with hole and 
electron pockets). More important for our purposes is the ac- 
tual value of the ordered moment for the iron spin, /zpe- Table 2 
of Ref. 8 and Table 2 of Ref. 9 provide a summary of the neu- 
tron scattering experimental values foryKpe: within the 1111 
and 122 families they range from 0.25 for NdOFeAs to ~ 1 
for SrFe2As2 and BaFe2As2 (the 1 1 family has substantially 
larger values of yUpe, but their ordering wavevectors are differ- 
ent, indicating that the 1 1 compounds require a special dis- 
cussion beyond the scope of this work). Since the mean-field 
investigations reported here focus on the mean-field solution 
corresponding to wavevector (tt, 0), the theoretical value of 
//Pe can be obtained when varying the couplings in the Hub- 
bard model, and a range compatible with experiments can be 
found. Not surprisingly, considering the (bad) metallic char- 
acter of these compounds, the range of interest will be shown 
to correspond to intermediate values of the on-site repulsion 
U, as also discussed before. ''^-^^ However, our investigations 
below show that sizable constraints over the Hund coupling J, 
not explored before, can also be obtained. 

Another powerful experimental technique that will help us 
to establish constraints on the couplings of the Hubbard mod- 
els is Angle Resolved Photoemission (ARPES). Applying this 
technique, several reports describing the Fermi surface of 
the pnictides have been presented, including their tempera- 
ture evolution. Here, the vast ARPES literature will not be 
comprehensively reviewed but the focus will be on the re- 
sults for the undoped compounds in the spin density wave 
(SDW) regime. In particular, it has been reported by several 
groups that in this magnetically ordered state the Fermi sur- 
face presents extra "features" in the vicinity of the F-point 
hole pockets that are not present in band-structure calcula- 
tions for the paramagnetic state. These extra features, mainly 
found for the 122 compounds, are described as the existence 
of "satellite" pockets caused by V-shaped bands (or by V- 
shaped band crossings). ^'^^ Our focus will be on trying to 
reproduce qualitatively with mean-field approximations these 
type of satellite extra features that ARPES investigations have 
systematically reported in the magnetic state. Note that near 
the original electron pockets at {n, 0) and (0, n) there are also 
several reports of interesting modifications of the paramag- 
netic band structure when at low temperatures. However, in 
our opinion, the ARPES analysis of the modifications to the 
band structure electron pockets appears less robust due to the 
presence of considerable noise in the raw data, more than for 
the case of the hole pockets. For this reason, the focus here 
will be on the existence of satellite pockets at low tempera- 
tures near the original F-point hole-pockets. Note also that 
these satellite features have been given an electron-like char- 
acter in some publications'*** but considering the uncertainties 
in ARPES studies, mainly caused by the considerable back- 
grounds, our analysis below will focus on finding regions in 
the U-J/U plane where any kind of extra pocket is induced by 
the magnetic order close to the original F-point hole-pockets, 
regardless of whether its character is electron- or hole-like. 



While finding constraints on U and J in multiorbital Hub- 
bard models is already interesting, our effort here continues 
with the analysis of the pairing states that are obtained in 
those "realistic" U-J/U regions, via the use of the Random 
Phase Approximation (RPA).^''^^ Exploring pairing tenden- 
cies is important since there is a growing controversy in the 
Fe pnictides investigations with regards to the symmetry of the 
superconducting order parameter. While a variety of photoe- 
mission experiments report the absence of nodesj^^ at least 
clearly with regards to the F-point hole pockets, other bulk 
experiments report results compatible with a nodal supercon- 
ducting state. "^^''° This controversy has not been resolved and 
it is one of the most important open problems in pnictides. 
Our investigations show that in the realistic U-J/U regimen 
identified here, for both three- and five-orbital models, the 
RPA dominant superconducting state is nodal (this includes 
the A ig state with zeros (or near zeros) of the order parameter). 
Moreover, within RPA there is a clear competition between a 
variety of states belonging to different irreducible represen- 
tations suggesting that different pnictides may have different 
symmetries with regards to the superconducting order param- 
eter All these competing states are nodal or quasi-nodal, at 
least within the limitations of our approximations. 

The organization of this manuscript is the following. First, 
the models and technique used are presented in Sec.|II] This is 
followed in Sec.|III]by the search for a realistic U-J/U regime 
for the three-orbital model. The case of five orbitals is pre- 
sented in Sec. |IV] The dominant RPA pairing tendencies in 
the realistic regime are presented in Sec. |V] Finally, conclu- 
sions are provided in SeclVll 



II. MODELS AND TECHNIQUES 

A. The model Hamiltonians 

In this effort, the three-orbital Hubbard model introduced in 
Ref. 3 1 will be used first. This model is purely based on the 
d electrons of Fe and it considers only the three orbitals d^^, 
dy,, and d„., widely believed to be the most relevant orbitals 
at the Fermi surface for the pnictides. The reader is referred 
to the original publication?' for a full description of the model 
and its band structure (which is in good agreement with ab- 
initio calculations). In momentum space, the model includes 
a tight-binding term defined as: 



/?TB(k) = J] 7^^"(k)<,,A,..: 
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2f2 COS kf + 2f 1 cos ky + 4f3 cos k^ cos A:,. - yu, (2) 

2f 1 cos k^ + 2t2 cos ky + Atj, cos k^ cos ky - fi, (3) 
2f5(cosfci- -H cosky) 

-\-Atf, cos kx cos ky - fi + A„, (4) 

r^' ^4t4 sin k, sin ky, (5) 

f^' = 2it^ sin kx + 4/f8 sin k^ cos ky, (6) 



TABLE I: Parameters for the tight-binding portion of the tliree- 
orbital model used here. The overall energy unit is electron volts. 
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0.02 0.06 0.03 -0.01 0.2 0.3 -0.2 -tijl 0.4 



J.23 _ jil _ 2if^ sjjj ^^ + 4,-fg sjjj ^^ j-Qs ^^ 
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where a bar on top of a matrix element denotes the complex 
conjugate, /i and v range from 1 to 3 and label the orbitals d_^, 
(1), dy; (2), and d„ (3). Since the Hamiltonian for a one-iron 
unit cell is here considered, then k runs within the correspond- 
ing extended Brillouin zone (BZ) -n < k^-, ky < n. The actual 
values for the hopping amplitudes are in Table I. 

The Coulombic interacting portion of the three-orbital 
Hamiltonian is given by: 

~ 27 ^ Si,(j, ■ Si,/j ^g-j 

+ J} (d. ^d. ,d.„,d.„.+h.c.), 

/_i^ i,Q-,T i,a,J, i,/3,J, i,AT ' 

i,a<fi 

where a,fi = 1,2,3 denote the orbitals, Si,Q. («i,») is the spin 
(electronic density) of orbital a at site i (this index labels sites 
of the square lattice defined by the irons), and the relation 
U' = U - 2J between the Kanamori parameters has been 
used.^' The first two terms give the energy cost of having two 
electrons located in the same orbital or in different orbitals, 
both at the same site, respectively. The second line contains 
the Hund's rule coupling that favors the ferromagnetic (FM) 
alignment of the spins in different orbitals at the same lattice 
site. The "pair-hopping" term is in the third line and its cou- 
pling is equal to J by symmetry. Note that the values used for 
U and J can be substantially smaller than the atomic ones, be- 
cause the interactions may be screened by bands not included 
in the Hamiltonian. The Coulombic interaction terms intro- 
duced above have been used and discussed in several previous 



publication; 



„13.18.20.25.31 



where more details can be found by 



the readers. All e nerg ies are provided here in electron volts. 
As shown in Ref. [3ll, the electronic density of relevance for 
this model is n = 4 to reproduce the expected Fermi surface 
in the paramagnetic regime. 

In the present investigation, two five-orbital models (also 
based only on the d electrons of Fe) have also been used, at 
electronic density n = 6. By supplementing the three-orbital 
model by more complicated five orbital versions, our main 
goal is to verify the self-consistency of our approach. In other 
words, if the many models, with similar Fermi surfaces by 
construction, would give quite different ranges of couplings 
for the compatibility with neutron and ARPES results, then 
this would raise concerns about the entire calculation. It turns 
out that, as shown below, the J/U and U ranges that are found 
to be physically reasonable are similar in all cases, demon- 
strating that our approach is self-consistent. With regards to 
the specific five-orbital models used here, the tight-binding 



parameters of one of them are in the Appendix, while another 
set of hoppings is from Ref. ,21, At t/ = 0, all these models 
provide a Fermi surface (see below) that compares well with 
experiments and band structure calculations for the 122 com- 
pounds. The Coulombic interactions for five-orbitals are the 
obvious generalization of the terms used for three-orbitals. 



B. The mean-field approximation 

To study the ground state properties of the models in- 
troduced before, a mean-field approximation will be ap- 
plied. This approximation was already presented in previous 
publications, '^-^^ but it is here also discussed for complete- 
ness. The simple standard assumption of considering only the 
mean-field values for the diagonal operators is followed;*^ 



('^i'^o-^j.v.o-') = \nfj + — cos(q ■ ri)mj6ijdfiv6o 



(9) 



where q is the ordering wavevector of the magnetic order. 
tifj and nifj are mean-field parameters (to be determined self- 
consistently) describing the charge density and magnetization 
of the orbital jj, respectively. The rest of the notation is stan- 
dard. Applying Eq. (|9]l to Him, the mean-field Hamiltonian in 
momentum space can be written as 



Hmf - H-YB + C + 2_j £iid^i^^d^,ii,a- 

+ 2 1n,<T(dli,,a''k+q,n,a + 4+q,/.,^^k,/.,<r), (10) 

ti,/i,a- 

where k runs over the extended first BZ, //tb is the hopping 
term in Eq. ([TJ, the constant C is 

C- - NuJ]lnl-^ml\-Ni2U'-J)Y,ft^fty 

11 ^ ' ll±V 
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A^ is the number of sites, and the following definitions were 
introduced 
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(11) 

(12) 



The mean-field Hamiltonian can be numerically solved for 
a fixed set of mean-field parameters using standard library 
subroutines. «^, and m^ are obtained self-consistently by min- 
imizing the energy via an iterative process. During the iter- 
ations 2^n;j=n was enforced at each step, such that the to- 
tal charge density is a constant (4 for the three-orbital model, 
and 6 for the five-orbital models). The reader should assume 
that these are the electronic densities used for these models 



throughout the manuscript, both in the mean-field approxima- 
tion and for the RPA approximation as well. Note also that 
the numerical solution of the mean-field Hamiltonian imme- 
diately allows for the calculation of the band structure, den- 
sity of states (DOS), and magnetization (m = 2/j "V) ^^ ^^^ 
ordering wavevector q. Moreover, the photoemission spectral 
function can also be calculated, as explained in Ref . I ISl 



ni. RESULTS FOR THE THREE-ORBITAL MODEL 

Our discussion of results starts with the three-orbital model. 
Some aspects of this discussion have been briefly mentioned 
in other publications, thus references to those previous efforts 
are provided where appropriate. 



A. Comparison with neutron scattering experiments 
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FIG. 1: Mean-field order parameter at wavevector (n, 0) vs. U (in eV 
units) for the three-orbital model discussed in the text, and parametric 
with the values of J/U indicated. 



Figure [T] contains the mean-field order parameter (m) at 
wavevector (n, 0) vs. U. The plotted values for m arise from 
the numerical solution of the mean-field equations discussed 
in the previous section. For small values of J/U, such as 0.00 
and 0.05, in Fig. [T]m discontinuously jumps from zero to a 
robust value at a critical U. While such a discontinuity is ob- 
served in all models discussed in this paper, its origin is not 
universal. It is caused by a metal-insulator transition in a four- 
band model, see Sec. II. C. 4 of Ref. J^. In the three-orbital 
model, the discontinuity only coincides with the opening of a 
gap for smaller J/U < 0.15, and is rather marked by the sud- 
den onset of strong orbital order, with a close competition be- 
tween substantial alternating and nearly perfect ferro-orbital 
order^' For larger a0.15 < J/U < 0.22, states with both types 
of orbital order can remain metallic for a small range of U 
just above the onset of strong orbital order, but the FS is qual- 
itatively very different from ARPES results, e.g. it does not 
feature any sign of hole pockets around the F point^i At or 
rapidly after the critical U, the density of states develops a 
gap (not shown), signaling insulating behavior, in contradic- 
tion with the experimentally observed (bad) metallic character 
of the undoped pnictides. Moreover, as a consequence of the 
discontinuity, the order parameter m in the range of small J/U 
never reaches the realistic values for pnictides reported in neu- 
tron scattering experiments, i.e. [0.25,1.0] for (n, 0) magnetic 
order Thus, these results for small J/U start illustrating one 
of the main messages of this publication, namely that within 
the mean-field approximation used in our effort the request of 
qualitative agreement with the experimental properties of the 
undoped pnictides imposes severe constraints on the values of 
U and J/U for the multi-orbital Hubbard models. In particu- 
lar, it is clear that J/U - 0.00 and 0.05 do not seem physically 
appropriate to describe the pnictides. 

As J/U increases further, m now develops (becomes 
nonzero) at an earlier critical value of U, allowing for a proper 
description of materials with weak magnetic order parameters 
such as the "1111" family. Note that for J/U - 0.10, 0.15, 
and 0.20 , a discontinuity is still present in m vs. U, so not all 
values of the order parameter m are possible, while for larger 



J /Us the m curves are no longer discontinuous. J/U - 0.50 is 
the largest ratio that should be considered to avoid a negative 
U' due to the relation U - U' + 27. Thus, adding this in- 
formation to the results for m and its comparison with neutron 
scattering, the proper range of J/ U couplings naively becomes 
~ [0.10, 0.50], with U larger than the first critical value where 
m develops. However, if in addition it is considered that J 
should be smaller than U' , then this reduces the range further 
to [0.10,0.33], since 7 = U'atJ/U^ 1/3. 



B. Comparison with ARPES experiments 

As discussed in the Introduction, another experimental 
source of information that can be used to reduce the allowed 
range of couplings in the Hubbard model is provided by the 
ARPES results for undoped pnictides. As previously men- 
tioned, a common generic feature of several ARPES experi- 
ments at low temperatures in the SDW phase is the develop- 
ment of "extra" features (pockets) near the original F-point 
hole pockets of the noninteracting limit. To search for these 
features, within our mean-field approximation the one-particle 
spectral function A(k, w) has been calculated, and the Fermi 
surface results have been analyzed in a wide range of U and 
J/U, using a ^-function broadening 0.025 eV. 

Shown in Fig. |2] are representative results of our ARPES 
calculations (see also Refs. 31,33). The focus is on a range of 
J/ U and U where mean-field ARPES contains a F-point hole 
pocket, as in the original C/ = bandstructure (presumably 
corresponding to the high temperature non-magnetic regime 
as well), and in addition "satellite" pockets as in ARPES, in 
between the original hole and electron pockets along the kjc 
axis. The results have not been folded, but are representative 
of a single-domain spin order wavevector, in this case (n, 0), 
and with only one Fe per unit cell. Figure |2] shows repre- 
sentative cases where these satellites are clearly present [the 
satellite pockets tend to be electron-like for small U, switch- 
ing to hole-like for slightly larger, but still realistic, values 
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FIG. 2: Unfolded mean-field Fermi surfaces for the three-orbital 
model corresponding to (m=order parameter) (a) (/ = 0, as reference 
(there are two hole pockets at F but they appear merged due to the 
broadening used for plotting); (b) J/ if = 0.33, U = 0.6, m = 0.2; 
(c) J/U = 0.20, U = 1.0, m = 0.4; and (d) J/U = 0.25, U = 1.05, 
m = 0.6. 



of U (for more details see Ref. QJ)]. The corresponding val- 
ues of J/U and U are indicated in the figure caption. These 
physically acceptable Fermi surfaces from the ARPES per- 
spective are found in the same approximate range of cou- 
plings as those selected from the m/neutrons perspective, as 
discussed above. It has been suggested before that a reduced 
ordered magnetic moment goes together with realistic A(k, w) 
in numerical approaches like density-functional theory^^ and 
the present mean-field scheme,^^ and we see here that this 
is not accidental, but that the analysis of both neutrons and 
ARPES are mutually consistent over a larger parameter range. 



as well as for a variety of models, see Sec.|lV] The four panels 
shown are qualitatively similar and further refinements in the 
so-called "physical region" (see Fig.lUi will need better tools 
for calculations and more accurate ARPES experiments. 
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FIG. 3: Unfolded band structure mean-field results for the three- 
orbital model and cases (a) C = 0, as reference; (b) J/U = 0.33, 
U = 0.6, m = 0.2; (c) J/U = 0.25, U = 1.05, m = 0.6. Panels 
(b) and (c) show a V-shaped pocket in between the (0, 0) and (n, 0) 
points. The scale used (arbitrary units) is on the right of the panels. 



In Fig. |3] some of the full spectral functions are shown, and 
compared with the U-0 case. The appearance of V-shaped 
features, that induce the presence of satellite pockets, is clear 
in these figures. These mean-field results for A(k, oj) are qual- 
itatively consistent with ARPES experiments for the pnictides 
that have reported similar V-shaped branchesi^ 

The results of Fig. |2] are in qualitative agreement with ex- 
periments, as already remarked in Ref. 22.- Moreover, our 
comprehensive analysis of A(k, a>) has shown that these fea- 
tures do not appear in other regions of the U-J/U phase di- 
agram. For instance, before the critical U where m develops 
from zero there are no satellite pockets, since they arise from 
the nonzero magnetic order and nesting effects. In the other 
extreme of U couplings larger than those used in Fig. |2] the 
Hubbard model simply becomes insulating (as discussed be- 
fore in Ref.USl). and there is no longer a Fermi surface. 



C. Summary phase diagram for three orbitals 

The values of the order parameter m and their compari- 
son with neutron scattering results, and the Fermi surfaces 
and their comparison with photoemission experiments, lead 
us to the mean-field phase diagram shown in Fig. IH one of 
the main results of this effort. In this figure, the range of U- 
J/U couplings compatible with neutron- ARPES experiments 
is labeled "physical region". This region is relatively small, 
providing substantial constraints on the parameters to be used 
in three-orbital Hubbard model investigations. If U is smaller 
than in the physical region, then the state is not magnetic; if 
U is larger, the state is insulating or it has a much distorted 
Fermi surface. If J/U is smaller than in the physical region, 
there is no room for the small and intermediate value order 
parameters found in neutron scattering; if J/U is larger, then 
U' becomes too small or negative and thus unphysical. 

It is important to clarify that our mean-field approximation 
does not incorporate the effect of fluctuations. For this reason 
this type of approximations are somewhat "rigid" and it can be 
expected,^^ although actual calculations are very difficult, that 
the true "physical region" may be larger than shown in Fig.|4l 
Thus, readers should consider the location of our "physical re- 
gions" as the center of a potentially broader area where agree- 
ment with experiments can be found. But this does not invali- 
date our main point: intermediate f/'s and intermediate J/U's 
are needed for agreement with available neutrons, transport, 
and photoemission experiments. 

A final remark is with regards to the actual value of U of 
order just 1 eV in the physical region: this cannot be the bare 
U but must already incorporate the influence of screening in 
a model where the long-range Coulomb interactions are in- 
cluded. It is for this reason that standard metals in general 
tend to have very small Us when studied via Hubbard-like 
models, while it is known that the bare atomic values for U 
are always of several eVs^ 



IV. RESULTS FOR FIVE-ORBITAL MODELS 

In this section, results for two five-orbital models are pre- 
sented, with a similar organization as for three orbitals. The 
hopping amplitudes and on-site energies of a novel "Model 
1" are provided in the Appendix (using as criterion for their 
determination finding qualitative agreement with band calcu- 
lations, as shown in Fig. 6(a)). The more accurate "Model 
2" is the model introduced in Ref. J2ll, where the reader can 
find the actual values of the parameters. These two models 
generate similar Fermi surfaces, but the values for the hop- 
ping amplitudes are rather different. Thus, they are useful to 
test whether our conclusions do or do not depend on small 
details. Indeed, an overall conclusion of our study is that the 
"physical region" is qualitatively similar for all the models an- 
alyzed in this manuscript. Both models studied in this section 
are at electronic density n - 6.0 in order to address the par- 
ent compoundsi^ The (5-function broadening used here is 0.01 
eV. 
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FIG. 4: Phase diagram for the three-orbital model obtained with the 
mean-field approximation described in the text. The "physical re- 
gion" in yellow is the regime of couplings found to be compatible 
with neutron and photoemission experiments. The "non-magnetic" 
region corresponds to a regime where the state has a zero order pa- 
rameter. In the "insulator" region, there is no Fermi surface and the 
state is insulating. The "discontinuity" label corresponds to the dis- 
continuous jump in the order parameter shown in Fig.[T] The entire 
"magnetic metallic" regime could in principle have been compatible 
with experiments, but only in the yellow highlighted region is that m 
is sufficiently small-intermediate in value and the Fermi surface has 
satellite pockets near the F-point hole pockets. 



A. Comparison with neutron scattering for Model 1 



The mean-field order parameter at wavevector (n, 0) vs. U 
for Model 1 is shown in Fig.|5] parametric with J/U. These 
results are in several aspects qualitatively similar to those dis- 
cussed before for three orbitals in Fig. [1] but there are some 
differences. For instance, in this case a discontinuity in m 
is found for all the values of J/U investigated. This shows 
that for this model and using mean-field techniques, there is a 
range of values of m (roughly between 1 and 2.5 depending on 
J/U) for which there are no solutions. It turns out that neutron 
scattering experiments for 1111 and 122 materials have un- 
veiled low to intermediate values of m, and such order param- 
eters fit in the range where the mean-field analysis provides 
stable solutions. Thus, for the case of five orbitals, the com- 
parison between the order parameter m and neutron scattering 
simply restricts U to be between the first critical value, where 
m develops, and the second critical coupling where the discon- 
tinuity occurs. This range is larger for small J than for larger 
values.^'' However, it should be noted that the magnetization 
in each orbital is parallel even for small J and U, in contrast 
to the state with antiparallel orbital magnetization found in 
Ref. I25I and in contrast to the model discussed in Sec. lIVE)] 
The small overall magnetic moment results here from weak 
magnetization of the individual orbitals rather than from their 
partial cancellation. 
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FIG. 5: Mean-field order parameter at wavevector (n, 0) vs. U (in eV 
units) for the five-orbital Model 1 discussed in the text, parametric 
with the values of J/U indicated. Note the presence of a critical 
U where the order parameter first develops, and a second critical U 
where a discontinuity occurs (see also Ref. flS) . 



B. Comparison with ARPES experiments for Model 1 

The comparison between the mean-field one-particle spec- 
tral function for the five-orbital models and ARPES exper- 
iments introduces more severe restrictions on the values of 
J/U and U than those discussed before in the neutron scat- 
tering context. Figure |6] shows the Fermi surfaces of Model 
1 for special values of J/U and f/ in a region compatible 
with ARPES. Qualitatively, these results resemble those of the 
three-orbital model: there are remnants of the original hole 
and electron pockets of the non-interacting limit and, in addi- 
tion, there are satellite features near the F-point (and also near 
the original non-interacting electron pocket at (n, 0)). 

Figure |7] shows the actual bands for the five-orbital Model 
1 for the non-interacting case and one example of a set of 
couplings with satellite pockets created by the magnetic or- 
der. While qualitatively similar to the results for three orbitals 
(Fig. [3]l there are some interesting differences: here the V- 
shaped features that originate the satellite pockets arise from a 
combination of two bands while for three orbitals they emerge 
from the bending of a single band (see also Ref.[33l). However, 
the level of accuracy of the ARPES experiments is not suffi- 
cient to distinguish between these two cases, plus there is still 
room to further refine the hopping amplitudes of the models 
used here to adjust for finer details of the ARPES experiments. 

Even in regions of parameter space where the ground state 
is magnetic and metallic, and the order parameter is in the 
range of neutron's experiments, the Fermi surface may still 
be qualitatively different from that observed in ARPES. As 
an example, consider the case shown in Fig. |8] correspond- 
ing to J/U - 0.10. These mean-field results for A(k, w) are 
clearly different from those shown before in Fig. |6] and re- 
gions where this type of discrepancies are found are removed 
from the "physical region" for the model. 
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FIG. 6: Unfolded Fermi surfaces of the five-orbital Model 1 in the 
mean-field approximation for the cases (a) 1/ = 7 = 0, as reference; 
(b) J/U = 0.23, U = 1.25, m = 0.2; (c) J/U = 0.28, U = 1.45, 
in = 0.7; and (d) J/U = 0.30, U = 1.4, m = 0.7. As in the case of 
three orbitals (Fig.|2]l and as in ARPES experiments, the results show 
distorted F-point hole pockets and satellite features next to them (and 
also next to the electron pockets at (n, 0)). 



C. Summary phase diagram for the five orbitals Model 1 

Similarly as for three orbitals, here a summary phase dia- 
gram is provided for the five-orbital Model 1 in Fig. |9] The 
labeling convention is the same as for three orbitals in Fig.|4l 
Using the information about neutron scattering and order pa- 
rameters restricts U and J/U simply to be between the non- 
magnetic region and the discontinuity line. From these per- 
spectives alone the "physical region" would be much larger 




FIG. 7: Unfolded band structure mean-field results for the five- 
orbital Model 1 and cases (a.) U = J = 0, as reference, and (b) 
J/U = 0.28, U = 1.45, m = 0.7, illustrating in (b) the origin of the 
satellite pockets shown in Fig. [6] The scale used (arbitrary units) is 
on the right. 
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FIG. 8: Fermi surface of the five-orbital Model 1 in the mean-field 
approximation for the case J/U = 0.10, U = 2.1, and m = 0.36. The 
lack of qualitative similarity with ARPES experiments shows that 
these couplings are not physically relevant to describe the pnictides. 
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FIG. 9: Phase diagram for the five-orbital Model 1 obtained with 
the mean-field approximation. As in Fig. |4] the "physical region" 
in yellow is the regime of couplings found to be compatible with 
neutron and photoemission experiments. The rest of the notation and 
details were already explained in Fig.|4] 



pared with the results for Model 1 in Fig. [Sj note that now 
the discontinuity in the order parameter occurs only at small 
J/U. In this respect, the results for Model 1 are similar to 
those recently reported for a similar model, ^^ where small J/ U 
and intermediate U were emphasized (see also Ref. 18). This 
similarity persists in the microscopic details of the state real- 
ized at small J: The values for the magnetization of different 
orbitals can here have a different sign, i.e., the intermediate- 
spin state with antiparallel orbital magnetization discussed in 
Ref.|25lis stabilized. While neutron scattering can only detect 
the overall moment and is not expected to distinguish between 
this scenario and the one with parallel magnetization found in 
other models, the precise microscopic nature of the state real- 
ized in iron pnictide compounds remains to be clarified. 



E. Comparison with ARPES experiments for Model 2 



than for the three orbitals case. However, considering the 
Fermi surface shape and its comparison with ARPES intro- 
duces more severe constraints, basically excluding the small 
7/ 1/ regime below 0. 15. As a consequence, the final "physical 
region" ends up being similar to that obtained with the three- 
orbital model. As explained in Sec. III.C, note also that fluc- 
tuations beyond our mean-field approximation are expected to 
render the "physical region" actually larger than displayed in 
Fig.ig 



D. Comparison with neutron scattering for Model 2 

The mean-field order parameter at wavevector (n, 0) vs. U 
is shown for Model 2 in Fig. [TO] parametric with J/U. Com- 



Similarly as presented for Model 1 and for the three-orbital 
model. Figure [TT] shows the Fermi surfaces for special val- 
ues of J/U and U in a region found to be qualitatively com- 
patible with ARPES. These Fermi surfaces show remnants of 
the non-interacting Fermi surfaces, plus satellite extra features 
near the F- and (n, 0))-points. Note that the character of these 
satellites, namely whether they are electron or hole pockets, 
depends on details and our focus has only been on the exis- 
tence of extra features in the correct location as compared with 
ARPES. Figure[T2]shows the bands for the five-orbital Model 
2, both for the non-interacting case and for one example from 
the set of couplings used in Fig. [TT] Most of the comments 
already made for the band structure of Model 1 with regards 
to the V-shaped features apply to Model 2 as well. However, 
it is interesting to notice that the satellite pockets are mainly 
electron-like for Model 1 and hole-like for Model 2. 




FIG. 10: Mean-field order parameter at wavevector (/r, 0) vs. U (in 
eV units) for the five-orbital Model 2, parametric with the values of 
J/U indicated. Note that compared with Fig.|5] the discontinuity in 
the order parameter now only occurs at small J/U. 



F. Summary phase diagram for the five orbitals Model 2 

As for the other models considered in this effort, in Fig.[T3]a 
summary phase diagram is provided for the five-orbital Model 
2. The labeling convention is the same as in Figs. |4]and|9l 
as well as the procedure to establish the so-called "physical 
region". The approximate location of this region with regards 
to U and J/U is similar to that reported in Figs. |4]and |9l 
showing again that our conclusions do not strongly depend on 
details. The physical region in Fig.[T3]is narrow along the U 
axis direction because m changes rapidly with increasing U, at 
small m. But, as discussed in Sec. III.C, fluctuations beyond 
mean-field approximations are expected to expand the size of 
the physical regions of the various models. 



V. RPA PAIRING SYMMETRIES ANALYSIS 

After having established a regime of couplings U and J/U 
where the mean-field approximation to the multiorbital Hub- 
bard model gives a qualitative agreement with neutron and 
ARPES experiments, and also with the metallic nature of the 
undoped compounds, it is interesting to investigate what kind 
of pairing tendencies are observed in those regions of param- 
eter space. Unfortunately, the number of many-body tools 
available to carry out such investigation is very small. While 
for the two-orbital model it is possible to carry out Lanc- 
zos studies on small clusters'-' to at least analyze the quan- 
tum numbers of the ground state with two extra particles, for 
three orbitals or more this calculation is no longer practical. 
There are also no sign-problem-free Monte Carlo techniques 
available for these complex calculations. Thus, here the Ran- 
dom Phase Approximation will be used. Experience with the 
cuprates indicates that this method did capture the relevance 
of (i-wave pairing in that context and for this reason it will be 
used here as well, although with the caveat that the method 
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FIG. 1 1 : Unfolded Fermi surfaces of the five-orbital Model 2 in the 
mean-field approximation at m = 6.0, for (a) U = J = 0, as reference; 
(b) J/U = 0.20, U = 1.35, m = 0.2; (c) J/U = 0.25, U = 1.35, 
m = 0.5; and (d) J/U = 0.30, U = 1.25, m = 0.6. As in the previous 
cases (Figs. 12] and |6j and as in ARPES experiments, the results show 
distorted F-point hole pockets and satellite features next to them (and 
also next to the electron pockets at (n, 0)). 



is qualitative at best. RPA has been recently applied to the 
five-orbital Model 2 considered here^' and the dominance of 
A[g states was unveiled at small J/U, with Big a close com- 
petitor In our investigations shown below, in the "realistic" 
regime of couplings the RPA dominant pairing tendency al- 
ways has nodes or quasi-nodes, an interesting result in view 
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FIG. 12: Unfolded band structure mean-field results for the five- 
orbital Model 2 and cases {a.) if = J = 0, as reference, and (b) 
J/U = 0.25, U = 1.35, m = 0.5, illustrating in (b) the satellite pock- 
ets shown in Fig.[TT| The scale used (arbitrary units) is on the right. 
Note the similarity of these results with those of Fig.|7] 



A. Random Phase Approximation Formalism 

For completeness, here the RPA approximation will be 
briefly reviewed, with emphasis on some technical aspects. 
This section follows the formalism already described in 
Refs.l2llandl65l Assuming that spin fluctuations (excitations 
in the paramagnetic state above the critical temperature, also 
called paramagnons) are responsible for the pairing mecha- 
nism present in the iron-pnictides, the RPA method will be 
used to study these fluctuations (caused by itinerant carriers) 
beyond the mean-field level. In linear response theory, prior 
to including the many-body interactions, these spin-wave-like 
excitations are obtained from the non-interacting Lindhard 
function (here defined for a multi-orbital model) 

xIimSI^ = ;^G,„,3(k + q)G,„,,(k), (13) 

k 

where each of the four indices needed to characterize the 
Lindhard function takes the values Z, = 1, ..., «„ (/ = 1, ...,4), 
with n„ the total number of orbitals being considered (here, 
rio - 3 or Ho - 5). The non-interacting Green's function, 
which describes the propagation of the elementary excitations 
present in the non-interacting model, can be written as 



G/,,/,(k,w) = ^ 



(^-EJk) 



(14) 
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FIG. 13: Phase diagram for the five-orbital Model 2 obtained with the 
mean-field approximation. As in Figs.|4]and|9l the "physical region" 
in yellow is the regime of couplings found to be compatible with 
neutron and photoemission experiments. The rest of the notation and 
details were already explained in Fig.|4] 



where b^lik) - (/i|//k) projects the tight-binding state lyuk) 
into the orbital |/i ), where // labels which one of the Irto tight- 
binding bands (in the extended BZ, with one iron per unit cell) 
the tight-binding state belongs to (with energy ^^(k)). In- 
troducing many-body interactions in the model, the Lindhard 
function gives origin to (distinct) spin and charge suscepti- 
bilities, which, calculated at the RPA level through a Dyson 
equation, can be expressed as (using matrix equations). 



Xs(q) 



Ml) 



r(q) 



1 - U'xHq) 



fig) 



(15) 



(16) 



where expressions for the non-zero matrix elements of U'' and 
{/^, in terms of the many-body interactions, can be found in 
Refs. 21 and ,65, The magnetic susceptibility (both RPA and 
bare) to be shown below in, e.g., Fig.[T4]is given by 



/T. 



:(q)=^2^.)(q);:;;;- 



(17) 



of the current experimental controversy between ARPES and 
other techniques with regards to the nodal structure of the su- 
perconducting state. It is also important to remark that other 
irreducible representations are closely competing with those 
that dominate, and small changes in parameters alter their rel- 
ative dominance. 



These RPA susceptibilities are used to construct a spin- 
singlet pairing interaction describing the exchange of charge 
(orbital) and spin fluctuations, resulting in an effective 
electron-electron interaction:— 

V\q) = lu'UqW' - ^U%(qW' + \(U' + U% (18) 
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Assuming that the dominant scattering occurs close to the 
Fermi surface, one can calculate the scattering amplitude 
of a Cooper pair between two points at the Fermi surface 
[(k, -k) — > (k', -k')], where the momenta k and k' are re- 
stricted to the Fermi surface pockets / and j, which span the 
existing pockets for the chemical potential used. To indicate 
this restriction a subscript / (or j) is added to the bands fi and 
V in question, and the (symmetrized) interaction vertex be- 
comes: 

ry(k,k') = Yj ^!r;*(-k)^!;;''(k)Refv;-;;(k,k')l 

h,h,h.U 
X b'^.(k')bl)(-k'), (19) 

which, after plugging into the linearized Eliashberg equation, 
results in a dimensionless pairing strength functional: ^^'^^ 



A[^(k)] 



1:0^,7^)4 Jog(k)r,.(k,kOg(kO 



(20) 



where the superconducting gap A(k) (in reciprocal space) was 
separated into an amplitude A and a function g(k), which 
should have the symmetry of one of the irreducible represen- 
tations of the corresponding point group (in the case of the 
pnictides, the D4/, group is considered, with representations 
Aig, A2g, Big, B2g, and Eg). The stationary condition leads to 
an eigenvalue problem defined over the Fermi surface: 



Z 



c/^"47r2vf(k') 



r/j(k,k')gj(k') ^ Agi(k), (21) 



which is an eigenvalue problem involving a matrix [r(k, k')], 
and where gi(k) is the value of the gap function for a point k 
at the pocket /, and A is the associated eigenvalue. The highest 
eigenvalue (normalized to 1) indicates what gap function will 
have the highest critical temperature. The dimension of the 
matrices to be diagonalized will be determined by the number 
of k points taken along the Fermi surface. For the calculations 
shown here, it was observed that considering approximately 
200 points along the Fermi surface provides a good conver- 
gence. All calculations were done at temperature T - 0.02eV, 
and an imaginary part 77 = 10^^ was used to regularize the 
Green's functions. All the sums over the BZ were done with 
uniform 64 x 64 meshes (calculations with a 128 x 128 mesh 
yielded qualitatively the same results). 



B. RPA Pairing Symmetries for the Five-Orbital Model 1 

The RPA analysis of the "realistic" U and J/U regime will 
start here with the five-orbital Model 1. In Fig. [14] the RPA 
results for the magnetic susceptibility xs are shown for two 
values of J/U (and slightly diff'erent U). The case J/U = 0.15 
is close but outside the "physical region" of Fig. |9] while 
J/U - 0.28 is clearly inside that region. However, with re- 
gards to xs, the figure shows that there are no substantial 
changes in the magnetic response for these two J /Us: an ap- 
proximately flat bare;^'5 becomes a RPA;^'^ with a sharp peak 
at the correct wavevector (n, 0). 
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FIG. 14: Bare and RPA magnetic susceptibilities vs. BZ wavevector, 
at the values of if and J/U indicated, for the five-orbital Model 1. 
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FIG. 15: (a) RPA pairing eigenvalues vs. U, at J/U = 0.15, for the 
five-orbital Model 1 . The symmetries of each eigenvalue are indi- 
cated, (b) Dominant A ig gap function with a similar color convention 
as in Ref. 21 (blue and red denoting different signs), (c) Subdominant 
state belonging to the Big representation. 



The RPA pairing eigenvalues A are shown in Fig. [15] (a) 
as a function of U, at a constant J/U = 0.15. Shown in (a) 
are the four first eigenvalues, that happened to have the four 
different symmetries that are indicated (as opposed to repeat- 
ing some symmetries). The dominant one is A [g, in agreement 
with Ref. 21. using the five-orbital Model 2. The actual pairing 
function for this dominant eigenvalue is shown in Fig.[T5](b). 
This is qualitatively of the form of the well-known s'^^ pair- 
ing, with different signs between the hole and electron pock- 
ets. However, the actual gap values are not uniform along 
each pocket. Actually, all the pockets present some naiTow 
regions where the gap functions nearly vanish, a feature quite 
reminiscent of nodes (dubbed here quasi-nodes*''). Overall, 
these results, and those presented below for Model 2 and for 
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the three-orbital model, are compatible with the recent obser- 
vations in Refs. I22ll65l that pockets at {n, n) suppress nodes in 
the A\g sector while the absence of those pockets produce a 
nodal A ig (note that in some of our results here and below, the 
superconducting order parameter at (tt, n) is weak and difficult 
to see). 

In Fig. [15] (c), the pairing function of the subdominant 
eigenvalue is shown. By mere inspection it is clear that it be- 
longs to the B\g sector, in qualitative agreement with Ref . l2ll 
that also identified this symmetry as the first competitor to A i^ . 
Thus, for JjU - 0.15 our results are similar to those previ- 
ously reported^' and have the value of providing a test of our 
procedure. However, Fig. [15] (a) contains extra information 
not discussed before: it also shows that two other eigenvalues 
are nearly degenerate with B\g and in view of the approxi- 
mate nature of the calculation they should also be considered 
as important competitors. These extra competitors belong to 
the Bog and Aog sectors. In particular the B2g is compatible 
with the Lanczos results found for the two-orbital model in a 
similar "intermediate U" region of parameter spaceJ^ 
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FIG. 16: (a) RPA pairing eigenvalues vs. U, at JjU = 0.28, for the 
five-orbital Model 1 . The symmetries of each eigenvalue are indi- 
cated, (b) Dominant Aig gap function with a similar color convention 
as in Ref. 21 (blue and red denoting different signs), (c) Subdominant 
state belonging to the Big representation. 



The presence of pairing states with nodes also extends to 
the regime of electron doping. For instance, for density n - 
6.125 our RPA studies (not shown) indicate the dominance of 
a B\g state for JjU - 0.28 and 0.35, with an Aig state with 
nodes being the subdominant pairing. While small changes 
alter the relative balance between the many competing states, 
all of those states, including the A \g, present a nodal (or quasi- 
nodal) structure at the RPA level used here. 

It is worth discussing these results in a historical context. 
In previous RPA studies of one-orbital Hubbard models in the 
spin-singlet channel, it was observed a dominance of the well- 
known iij.2_,,2-wave state over other channels.*^ Calculations 
in the one-orbital context did not report the presence of other 
pairing states so close to the dominant one as shown here and 
in Ref. j2ll. In fact, the only competitors to the li-wave state 
that have been identified for the case of one orbital are spin- 
triplet i-wave odd-frequency states (not studied in our present 
investigations due to their spin-triplet nature). The appear- 
ance of odd-frequency states as competitors was analyzed in 
studies on square lattices,^" as well as on triangular lattices. ^i 
But among the spin-singlet states, the dominance of d-wave 
was clear for one orbital. Thus, finding in the multiorbital 
models so many spin-singlet even-frequency states compet- 
ing with the dominant ones is surprising and merits further 
work. These results suggest that within the pnictide family 
different pairing channels may be stable in different pnictide 
compounds since small changes in parameters, such as caused 
by chemical doping or pressure, may lead to changes in the 
dominant pairing tendency. 
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FIG. 17: Bare and RPA magnetic susceptibilities vs. BZ wavevector, 
at the values of U and 7/(7 indicated, for the five-orbital Model 2. 



Figure [T6l provides novel interesting information. Switch- 
ing now to J/U - 0.28, i.e. inside the U-J/U region con- 
sidered realistic based on neutron and ARPES experiments, 
the dominant eigenvalue now belongs to the A2g irreducible 
representation (with a gap function shown in (b)). Big (see 
(c)) is still the subdominant tendency, but it is nearly degen- 
erate with Z?2g. For this value of J/U the Aig channel is 4th 
in the relative order Thus, a relatively small variation in J/U 
induces a qualitatively drastic reaiTangement of pairing eigen- 
values. However, all of them have nodes (or quasi-nodes as in 
the case of Aig). Thus, it appears unavoidable to conclude that 
nodal superconductivity should be present in the pnictides, at 
least within pairing tendencies that rely on electronic mecha- 
nisms (and within the RPA and mean-field approximations). 



C. RPA Pairing Symmetries for the Five-Orbital Model 2 

For completeness, the RPA analysis for Model 2 is here 
included. Model 2 was already studied in Ref. |2l], but here 
the focus will be on the "realistic" regime of U and J/U. In 
Fig. [17] the RPA results for the magnetic susceptibility ;t^5 are 
shown for representative values of U and J/U. In this regime, 
there is a clearly dominant peak at wavevector (n, 0), although 
there are subdominant peaks at other momenta as well. The 
results at other value of U and J/U in the vicinity of the one 
shown are similar. 

The corresponding RPA pairing eigenvalues A are shown in 
Fig. [13(a) as a function of U, at a constant J/U - 0.28. Fol- 
lowing a similar organization as in the case of Model 1 , shown 
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FIG. 18: (a) RPA pairing eigenvalues vs. U, at J/U = 0.28, for the 
five-orbital Model 2. The symmetries of each eigenvalue are indi- 
cated, (b) Dominant A i^ gap function with a similar color convention 
as in Ref. 21 (blue and red denoting different signs), (c) Subdominant 
state belonging to the Big representation. 



in (a) are the four first eigenvalues. The dominant one is Aig 
with nodes, in agreement with Refs. l2lll22L The pairing func- 
tion for this dominant eigenvalue is shown in Fig.[T8](b): it has 
different signs between the hole and electron pockets, and the 
actual gap values are not uniform along each pocket, actually 
presenting nodes. In Fig. [18] (c), the pairing function of the 
subdominant eigenvalue is shown. It belongs to the Big sec- 
tor, as found in Ref. 21. The importance of our results is that 
using exactly the same formalism our approach does properly 
reproduce previous results with regards to the dominance of 
the Aig sector for Model 2, while Model 1 (with a very similar 
Fermi surface) has pairing in other channels. This highlights 
again that very small changes in parameters can dramatically 
alter the pairing instability channel. 



D. RPA Pairing Symmetries for the Three-Orbital Model 

To complete the RPA analysis, now the three-orbital model 
will be considered. The results reported below will not be as 
clear as in the previous case of five orbitals, but these results 
are presented here anyway to alert the reader of the subtleties 
associated with RPA approximations. In spite of the difficul- 
ties to be shown below, the pairing states that dominate still 
have nodes. Thus, the general conclusion that nodal super- 
conductivity tends to be favored in Hubbard models remains 
the same, at least within the approximations used in our effort. 

Figure [19] shows the RPA magnetic susceptibility for two 
cases of interest, J/U - 0.15 and 0.30, within the "phys- 
ical region" of Fig. [4] These results already illustrate the 
main problem found here: although there is a small peak at 
wavevector (n, 0), there are other larger peaks at wavevectors 
closer to (0, 0). Thus, the results to be shown for the pairing 
tendencies are in a magnetic background that does not cor- 
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FIG. 19: Bare and RPA magnetic susceptibility vs. BZ wavevector, 
at the values of U and J/U indicated, for the three-orbital model. 
Note that in addition to a small peak at (n, 0), there are larger peaks 
at other wavevectors. 
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FIG. 20: (a) RPA pairing eigenvalues vs. U, at J/U = 0.25, for 
the three-orbital model. The symmetries of each eigenvalue are indi- 
cated, (b) Dominant B2g gap function with a similar color convention 
as in Ref.|2l| (blue and red denoting different signs), (c) Subdominant 
state belonging to the Aig representation. 



respond to that of the pnictides experiments, but it contains 
states with a variety of wavevectors. 

Figure [20] (a) displays the U dependence of the pairing 
eigenvalues at J/U = 0.25, inside the "physical region" for 
three orbitals. Note that in this case it is the B2g symmetry that 
dominates, as in Ref.[l3l with Aig being subdominant. The or- 
der parameter of the dominant B2g tendency is in Fig.|20l(b), 
while the subdominant A i^ is in Fig.|20](c) (with the exotic de- 
tail that in this case the two hole pockets caiTy a different sign 
for the order parameter). The presence of the B2g tendency 
was not only observed at the couplings of Fig. [20] and n - A, 
but also in a wide range of J/U and varying the electronic 
density n to 4. 1 and 3.9 (not shown). Thus, such nodal pairing 
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tendencies are robust. However, a more detailed analysis of 
the RPA response of the three-orbital model, motivated by the 
multiple wavevectors in its magnetic state, should be carried 
out in the future. 
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VI. CONCLUSIONS 

In this publication, the undoped limit of the multiorbital 
Hubbard model for pnictides has been studied using a stan- 
dard mean-field approximation, similar to that employed in 
the study of the single-orbital Hubbard model for the undoped 
cuprates. Within this approximation, the magnitude of the or- 
der parameter associated with the (n, 0) magnetic order was 
studied varying U and J/U. In addition, the one-particle spec- 
tral function A(k, co) was also analyzed. Comparing results 
against neutron scattering and ARPES experiments, allow us 
to define regions in parameter space, dubbed "physical re- 
gions", where the mean-field model Hamiltonian predictions 
are in qualitative agreement with the above mentioned exper- 
iments. These regions are relatively small in size since the 
ground state in this regime must be simultaneously metallic, 
magnetic with order parameters in the range found by neu- 
trons, and with Fermi surfaces containing satellite pockets in- 
duced by the magnetic state near the F-point hole pockets of 
the original band structure. Although fluctuations beyond the 
mean-field approximation are expected to enhance the physi- 
cal regions, our results still provide important constraints on 
the couplings to be used for theoretical studies of multiorbital 
Hubbard models for pnictides. 

In addition, in the regime of U and J/U described above, 
the RPA approximation allowed us to make predictions about 
the dominant pairing tendencies in the "physical regions". 
While it is clear that several channels are competing in these 
models, namely that small changes in parameters can lead 
to drastic changes in the dominant symmetries of the pair- 
ing states, the common property that emerges is the presence 
of nodal superconducting states in those physical regions. 
Within these models it appears very difficult to stabilize states 
without nodes. Even in regimes where the Aig state domi- 
nates, it still has quasi-nodes (very small values of the ampli- 
tude at particular Fermi surface points) or true nodes. Thus, 
in these regards our results are more compatible with the bulk 
measurements that reported nodal superconductivity than with 
the ARPES experiments reporting nodeless superconductivity 
in doped pnictides. However, our studies are based on approx- 
imations that need to be refined. Hopefully, our study will ini- 
tiate a debate on what are the true dominant tendencies in the 
multiorbital Hubbard model for pnictides, helping to decide if 
the mechanism is electronic or phononic. 
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Appendix: parameters of five-orbital Model 1 

Using a similar notation as for three orbitals, and the pa- 
rameters of Table II below, the tight-binding portion of the 
five-orbital Model 1 is HTB(k) = 2k,<7-,/j,v ?'/i,v(k)t/[_^_^t/^ „^ , 
where 

Tl 1/22 = ^tl^/y COS kjt + 2tyl^ COS ky + 4tly cos k^ COS ky 

±2tll(cos 2k^ - cos 2fey) + 4-t\l^,^^^,^, cos 2fct cos ky 
+4^ "y/xxy COS kj, COS 2ky + 4f "^,^, cos 2kj, cos 2ky 
+eii/22, (22) 

r33 = 2tl^ (cos kx + COS ky) + 4-tll COS k^ COS ky 

+2tll(cos 2k, + cos 2yty) + 633, (23) 

T44 - 24'*(cos kx + cos ky) + 4-t'l* cos k, cos ky 



+2t^x{cos 2k, + cos 2ky) 



-i-4f^^ (cos 2kj: cos ky + cos k^ cos 2ky) 



+4tjj^yy COS 2kj: COS 2ky + 644, 



(24) 



T^s - 2f^^(cos kx + cos ky) + 2fJ^(cos 2kx + cos 2ky) 
-i-4f (cos 2kj! cos ky + cos k^ cos 2ky) 

+4tfxyy COS 2kx COS 2ky + 655 , (25) 

7^12 = 7^21 



- -44y sin kx sin ky 

—4-tl^y(sm 2kx sin ky + sin k, sin 2ky) 
-4/ , sin 2kx sin 2ky, 



(26) 



7^13/23 - 7^31/32 



±2it^x sin fcy/A- ± 4/f J.y sin fey/^ cos k^iy 



+4/fj,^^,(sin 2kyjx cos k^jy - cos 2kxjy sin kyj^, (27) 

ri4/24 = Ta\JA1 

- 2it]* sin kx/y + 4itly sin kx/y cos ky/x 

+Ait\^xy sin '^kxiy cos kyix, (28) 

ri5/25 = 7^51/52 

= 2it]^ sin kyix - 4-itly sin ky/x cos kx/y 

-4it\\yy sin 2ky/x cos 2kx/y , (29) 

r34 = r43 

- 4f^^y(sin kx sin 2ky - sin 2kx sin ky), (30) 

7^35 = 7^53 

= 2f^ (cos kx - cos ky) 

-i-4f^^ (cos 2kx cos ky - cos kx cos 2ky), (31) 

Tm - TsA 

- 4txy sin kx sin ky + 4txlyy sin 2kx sin 2ky. (32) 
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TABLE II: Parameters for the tight-binding portion of the five-orbital 
Model 1 used here. The overall energy unit is electron volts. 
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